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Abstract: Recent SPS data on the rapidity distribution of protons in p+S, p+Au and 
S+S collisions at 200 AGeV and preliminary Pb+Pb collisions at 160 AGeV are compared 
to HIJING and VENUS calculations as well as to predictions based on the Multi-Chain 
Model (MCM). The preliminary Pb data suggest that a linear dependence of the proton 
rapidity shift as a function of the nuclear thickness, as first observed in p+A reactions, may 
apply up to Pb+Pb reactions. The observed rapidity dependence of produced hyperons 
in both p+A and A+A reactions however cannot be explained in terms of such models 
without introducing additional non-linear effects. 


1 Introduction 


Preliminary data on baryon number transport in Pb + Pb —> p,p,A,A i 1- I, I, i at 
160 AGeV has become available. These data are of interest as tests of nuclear stopping 
power ||. The first estimates of the baryon stopping power of heavy nuclei was made by 
Busza and Goldhaber [^] based on the A dependence of p + A —»■ p + X data at fixed 
Pj_ = 0.3 GeV. Baryon stopping power refers to the transport of baryon number away 
from the nuclear fragmentation regions and is measured in terms of the single inclusive 
rapidity distribution of protons and hyperons. In ref. [|J, the estimates were refined by 
taking into account exact Glauber geometry and applying the Multi-Chain Model (MCM) 
parameterization of baryon number transport. For a review of the average rapidity loss 
for p + p, and A + A collisions at beam momenta 11.6, 14.6 and 200 GeV/c per nucleon 
see ref. |Q. Until recently, the most complete information on nuclear stopping power was 
limited to the p + Ag data of Toothacker et al. 
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Systematic data from heavy-ion collisions [I], []. [Tj. Q [II], [f^. [I3j| provide new informa¬ 

tion about the nuclear stopping power. Important information about baryon inelasticity 
also comes from the analysis of veto and transverse energies distributions [|I4|| -|1I8| . 


As summarized in Q, baryon number transport is one of the key observable that has 
been debated for some time in connection with the ongoing search for nonlinear dynamical 
phenomena in nuclear reactions. One source of nonlinear behavior may arise if a quark- 
gluon plasma is formed in such reactions. In this connection, the observed strangeness 
enhancement and hyperon production data [[I] Q [II], [L4|,[|nj-[^2] have stimulated par¬ 
ticular interest |^3|, Q. In |25|] it was pointed out, however, that the unusual difference 
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between hyperon production in p + p and minimum bias p + S may reflect more the onset 
of novel non-equilibrium dynamical mechanisms. 

The breakdown of the linear dependence of baryon rapidity shifts as a function of the 
mean collision number, v A 1 / 3 , is one of the obvious places to look for new phenomena. 
The linearity assumption underlying MCM and many Monte Carlo event generators is 
motivated by the phenomenology of string formation and fragmentation PBL [ 




Such linear behavior is also expected if the radiative energy loss per unit length is a 
constant in nuclei . 

Several recent theoretical developments have raised, however, new possibilities for 
nonlinear dynamics associated with baryon number transport. In one development |31 


a non-linear energy loss of gluons in nuclei was predicted as a function of the nuclear 
thickness, L 

A E, « ±a,N e L (pp t ~ 15 GeV (AC) . (1) 

This energy loss could, under appropriate conditions, lead to enhanced rapidity shifts of 
hadrons passing through nuclei. 

A second possible source of nonlinearity was suggested in ref. [[52] based on the Rossi- 
Veneziano ||33|| baryon-junction Regge exchange model. This mechanism can give rise 
to large rapidity shifts of the baryon number offset by only a modest enhancement of 
the forward final pion rapidity density. This mechanism could even modify significantly 
baryon transport up to RHIC energies because the assumed junction Regge intercept, 
«j(0) = 0.5, leads to a slow suppression of this mechanism with increasing energy (oc 
s“ 1/4 ). 

An even more exotic baryon transport mechanism was proposed in [E3J based on a 
variant |35| of the above baryon-junction exchange model. In that variant the junction 
trajectory is assumed to have unit intercept and thus leads to an energy independent 
uniform rapidity density component of the inclusive proton yield. 

Given the new data and the above theoretical speculations, it is appropriate to take 
analyze carefully the available information on baryon transport and hyperon production. 
We therefore take into account the older p,p + Ag data at 100A GeV [[jUj, the new p + A 
0 at 200 GeV, the S + A at 200 AGeVQ, 0 [H 

i i !, 
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and the Pb + Pb data at 160A 
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In a previous paper, we concentrated on the anomalous hyperon 
Here we concentrate on the proton rapidity shifts in nuclear 


GeV 

production processes 

collisions. We test whether nuclear stopping power extrapolates linearly from p + A to 
Pb + Pb. In this analysis we recall the predictions of the Multi-Chain Model|8] (MCM) 
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as well as utilize the Monte Carlo HIJING1.3 


and VENUS4.12 [291 models as in 12 


2 The Multi Chain Model 


Since the HIJING and VENUS models were discussed extensively in |25|], we recall here 


only the essential elements of the MCM to be used in the present analysis. 

The single inclusive proton rapidity distribution in B + A —> p + X reactions is given 
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in general by Glauber theory as 

]N BA^pX B A 

——7-= 51 55 p BA{m,n)Q m , n (y) , (2) 

a y 171=1 n=1 

where the tedious but well understood nuclear geometry is separated from the multiple 
collision dynamics encoded in the functions, Q m , n (y)- The probability that a group of m 
projectile nucleons multiple scatter with n target nucleons is given by 

/' r/ 2 b /* d 2 s 

P B A(m,n) = /- B(b) / — P B (m,N B (s)/B)P A (n,N A (b-s)/A) , (3) 

where P A (m,x ) = C m:A x m ( 1 — x) A ~ m is the binomial distribution. The mean number 
of inelastic collision in nucleus A at impact parameter b in terms of the diffuse nuclear 
density and inelastic pp cross section, cr in ~ 32 mb by N A (b) = a in f dz p A (z, b). The 
impact parameter prohle function B is included above to account for the experimental 
trigger bias (e.g. veto or transverse energy cuts). For minimum bias events, B = 1. For 
a given prohle, the total reaction cross section is 

= J d 2 bB(b) (l - (1 - / d 2 sN B (s)N A (b - s)/(ABa in )) BA ) . (4) 

The main simplifying assumption of the MCM is that of independent fragmentation: 

Qm,n(y ) ~ ™Qn{Y ~ y) + nQ m {y) , (5) 

where Y is the rapidity difference between the projectile and target. As emphasized in ||, 
this linear superposition ansatz is far from obvious, and one of the important questions 
awaiting the recent heavy ion data is whether this breaks down for sufficiently heavy 
nuclear collisions [6[|. With this ansatz the rapidity distribution simplifies to 

dN BA ^P x A B 

——7-= r B W B 55 PB/A(n)Q n (Y - y) + r A W A 55 PA/B(m)Q n (m ) , (6) 

dy n =1 m=l 

where r B = ( Z B f + N B (1 — f))/B is the fraction of projectile baryons that fragment into 
protons, and where W B is the average number of wounded baryons in nucleus B :, 

W B = f —B( b) [ —N b ( s){1 - [1 - N a (b - s)/A] a } (7) 

J OAB J ®in 

and P B / A (n ) is the fraction of them that interact with n target nucleons. Typically 0, 
the fraction of incident protons that remain protons is / ~ 0.53 away from the diffractive 
peaks. 

What remains then to be determined are the dynamical fragmentation functions, 
Q n (y), specifying the rapidity distribution of target baryons that have suffered n inelastic 
interactions. For n — 1, i.e., p + p, the MCM assumes the simple form 

Qi{y,<x) = e~ y , ( 8 ) 
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as consistent with the flat dN/dx distribution observed away from the diffractive peak 
(x ~ 1). For n > 2 the multiple collision contribution to baryon rapidity transport is 
parameterized in the MCM model by a one parameter class of functions: 


Qn(y,a) = 


a 


a 


n —i 


n —2 


?-y 


D -ay 


£ —(a -1 )”y 


m =0 


ml 


( 9 ) 


These functions arise using a scaling algorithm|| in which the probability density that 
spectator partons retain a fraction z of the total light-cone momentum after an inelastic 
interaction is az “ -1 . This ansatz implies a geometric scaling of the spectator energy 
fraction moments (z p )Fn = («/(a + p)) n and a linear scaling of the mean rapidity shift 
with collision number n: 

( y)n = 1 + {n- 1 )/a . (10) 

Detailed fits to the Barton et al|! p + A —> p + X data at 100 GeV, gave a — 3 ± 1. 
Other models can also achieve good fits to the existing p + A data with different pararn- 
eterizations. For example, a recent fit with the constituent quark model was presented 
in [[I?]]. The advantage of the MCM approach is in its simplicity, reducing the problem 
of the nuclear stopping power to one phenomenological parameter, a. Recall that the 
naive incoherent cascade limit corresponds to a = 1, leading to one unit rapidity shift 
per interaction. The empirical a = 3 arises because the inelasticity per interaction inside 
nuclear targets is reduced by the finite formation time of secondary fragments. 

We note that the scaling approximation used above certainly breaks down for lower 
energies where the rapidity gap, Y, between projectile and target becomes ~ 3. In the 
following we simply cutoff Q n with a 9(Y — y ) factor and normalize Q n in that finite 
interval. In [§] a slightly different finite energy cutoff prescription was used. Fortunately, 
for Ei a b > 100 GeV such cutoff effects are unimportant. 

At the pp level we can test for the possibility of novel baryon junction transport |32[ . 
If the probability of junction exchange is ej and its Regge intercept is aj, then Q 1 would 
be modified from (BD to 


Qi(y) = (1 - ej)e y + ej(l - aj)e (1 aj)y 


( 11 ) 


Actually, the contribution of junction exchange close to the fragmentation regions is not 
well determined and the above form may only apply for ?/ > 1. Given the above form, 
the final p — p rapidity distribution in the cm with a rapidity gap 2 Y would be 


dN pp ^ pX 

dy 


= t p (1 - e j) 


cosh (y) 
sinh(V) 


+ ej( 1 - aj) 


cosh((l - aj)y) \ 
sinh((l - otj)Y)) 


( 12 ) 


In figure 1, the valence proton, p — p, distribution from pp collisions at 400 GeV[38 


is 


shown compared to MCM and HIJING. These same data were used to test VENUS in |29| . 
The mid-rapidity yield is underestimated by HIJING by a factor of two, while the VENUS 
proton fragmentation scheme fits better j29| . Both VENUS and HIJING over-predict the 
forward y ~ 6 yields. The excellent agreement with the simple MCM form, cosh (y — Y), 
with r p = 0.53 leaves little room for exotic baryon exchange contributions at this energy. 
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At higher energies, y/s = 53 AGeV, Ref. [[3^] finds evidence for taking ej > 0. However, 
Fig. 1 demonstrated conventional baryon trajectory exchange reproduces accurately the 
observed valence proton rapidity distribution at this lower energy. This point is important 
since in |J2| it was hinted that at SPS energies such conventional mechanism at least as 
expressed in the Dual Parton Model [27] could not reproduce the experimental stopping 


power in nucleus-nucleus reactions. Part of the problem in DPM and HIJING discrepancy 
with pp —> pX data can be traced to the assumed diquark fragmentation schemes. We 
return to this point in section 5. 


3 Baryon transport in p + A 

In Figure 2. we compare calculations with the 100 GeV p,p + Ag —> pX dataffTOfl . Part (a) 
shows the dN/dy distribution in the projectile fragmentation domain (A y < 2) and the 
higher but narrower dN/dy distribution in the target fragmentation domain (A y > 4). 
Both the MCM prediction (with a = 3) and the VENUS model are seen to reproduce 
the leading proton distribution within the experimental errors. HIJING leads to a sim¬ 
ilar mean rapidity shift but is distributed more narrowly about the mean. This is due 
to the default diquark fragmentation scheme in JETSET|26| which is used in HIJING. 
VENUS does not use the JETSET scheme and has been adjusted to reproduce the flat 


fragmentation region ||29 


The target region is isolated in part (b) through the p + A —> p + X channel. Since 
MCM and HIJING only account for the wounded baryons, the narrow nuclear enhance¬ 
ment around the target rapidity (A y = 5.3) can only be reproduced by VENUS, which 
incorporates final state interactions. However, the main interest is the probability that 
a target proton emerges near the projectile rapidity Ay ~ 1 in Fig.2b. While the data 
fluctuate greatly in that region, they are consistent with the expected exponential fall 
from eq.(^). As in Figure 1, there is little need here to invoke a junction exchange 
contribution |32| with intercept 1/2 that would lead to a exp[y/2] rather than the conven¬ 
tional exp[y] tail in this opposite fragmentation region. 

In Fig. 2c the projectile fragmentation region is isolated through the p + A —> p + X 
channel. The MCM model is consistent with the data while both HIJING and VENUS 
tend to overestimate this distribution just as both over-predicted the y ~ 6 region. In Fig. 
2d the A dependence of the MCM baryon distribution is illustrated. We recall that in 
MCM, dN/dy(y —► 0) = v p Pa{ 1) is completely fixed by the geometrical probability that 
the projectile proton interacts only once and the probability, r p that the proton remains 
a proton after fragmentation. (As in Fig.l , r p = 0.53 || here.) 

In Figure 3, we compare the models to preliminary p + A data reported from NA35|I|. 
Parts (a) and (c) show the p — p distributions for S and Au targets respectively. The 
peculiar feature in (a) relative to previous data is that in this case HIJING best reproduces 
the projectile fragmentation peak in p + S. In this case the MCM distribution is too flat. 
The mean rapidity shift in all three models is about the same (A y sa 1.3) in accord with 
previous systematics [0, but the narrow peak at y = 5 is unexpected for minimum bias 
events given in Figures 1 and 2. The new p + Au minimum bias data in part (d) also 
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indicates a greater stopping power than encoded in the models HIJING and VENUS with 
a significantly suppressed yield beyond y > 4 and an enhanced yield below y < 4. The 
calculated minimum bias rapidity shift in MCM is Ay ps 2 while the data indicate perhaps 
Ay « 2.5. Such a large rapidity shift is achieved only in the extrapolated most central 
p + A reactions in the older data [fl6H . 

There appears therefore to be a discrepancy between the data in Fig.3a and the older 
data|36|. Unfortunately the data in Fig.2 have too low statistical significance to settle 
this problem. Note that in these plots the unobserved target fragmentation regions are 
artificially cutoff below y < 0.5. Upcoming p + A data with NA49 will hopefully clarify 
this ambiguous situation. 

As seen in Fig. 3, the net A — A hyperon distributions is even in greater disagreement 
with respect to both VENUS and HIJING calculations. In ref. f25| the absence of A 


fragments beyond y > 5 in p + S and beyond y > 4 in p + Au was emphasized to be 
anomalous relative to pp —> A data and pp —> A data previously analyzed in [^5], |9f . In 
those p target reactions, the A rapidity distribution closely mirrors the proton fragment 
distribution with reduced normalization due to the expected suppression of strange flavor 
production. The suppression of the forward A relative to forward p cannot be accounted 
for even with the double string mechanism added into VENUS. The lack of any measured 
A’s beyond y > 4 in minimum bias p + Au is especially peculiar. We thus re-confirm that 
hyperon transport in p + A is anomalous relative to non-strange baryons, and is already 
evident in minimum bias p + S reactions, where on the average, the incident protons only 
interacts with two target nucleons. 


4 Baryon Number Flow in A + A 

We turn next to the distribution of valence baryons in A + A collisions at SPS for which 
some preliminary data have become available from NA49[^. [I|. |T2|] and NA44[J3J. As noted 
in || data on nuclear stopping power in Pb + Pb collisions is long awaited as a critical 
test of nuclear transport models and to constrain the maximal baryon densities achieved 
in such reactions. 

Figure 4 compares the spectrum of pion and participants protons in S + S at 200 

m 

7T~ 

(see reference 


AGeV 


The various data sets from 


and Pb + Pb reactions at 160 AGeV []2], [| jl 
NA35 for S + S —> tt~ + X correspond to different centrality triggers [|20|, with the higher 
one corresponding to a more severe veto trigger cut (see reference 


for more details). 

We see that the negative pion rapidity densities are well accounted for by both HIJING 
m and VENUS 


291. This is largely due to the fact that the pion distribution simply 


grows linearly with the atomic number between S + S and Pb + Pb. The centrality trigger 
was implemented in the above calculations via an impact parameter cut of b < 1 fm. For 
the MCM calculation the impact parameter cut for Pb + Pb was taken to be 3.3 fm to 
simulate more closely the 5-6% centrality trigger |2|. 

The significant central rapidity minimum in the p — p dN/dy predicted by HIJING is 
traced back here to the excess central minimum in pp in Fig.l characteristic of FRITIOF 
type string models. By adjusting the fragmentation functions that fit the pp central 
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rapidity region better, VENUS, can avoid the suppression of the mid-rapidity protons 
that HIJING and other models using JETSET predict. On the other hand, MCM, which 
fits well both the pp and the pA data as shown in Figs. 1,2, is seen to reproduce well the 
central p — p distribution in both SS and PbPb. 

Figures 4a,c demonstrate the insensitivity of the pion distribution to the underlying 
baryon number flow. In particular, the rather large difference between the proton distri¬ 
bution in HIJING and VENUS in Fig.4d is contrasted by the much more modest difference 
of the pion distributions in Fig.4c. As we show in section 5, the forward energy flux is a 
more sensitive measure of the energy degradation difference between the models. 

Our main conclusion is that the central region baryon number transport in A + A 
can be well understood on the basis of a linear extrapolation of the mean rapidity shift 
as a function of the collision number as given by eq.(fl|). The main difference between 
the valence distribution functions in pA and BA arises simply from the variation of the 
Glauber geometrical probabilities of multiple interactions in the different projectile and 
target combinations. 

An important caveat to the above conclusion is the somewhat narrower rapidity dis¬ 
tribution of the most recent analysis Q on Pb + Pb (p — p) than reported in [|2J. This 
difference is shown by comparing the solid dots to the solid triangles in Fig.4d.If in the 
final analysis, the [Q narrow distribution is conhrmed, then the Pb + Pb baryon stopping 
power would have to be regarded as anomalous relative to p + A|7], |T(|has well since then 
an additional 1/2 unit rapidity shift would be required over that predicted by MCM nor 
VENUS. In [|4[, it was estimated that the net integrated proton number implied by the 
recent analysis is close to 2 x Z(Pb). This itself is peculiar since the one expects more 
protons from n —> p exchange. At this time the central rapidity region is most reliable 
since another experiment [[3] and alternate TOF measurements converge within errors in 
that region and are closest to the theoretical expectations based on p + A systematics. 
Nevertheless, it will be important to follow up the high rapidity tails assess the significance 
of the present preliminary indications. 

Figure 5 shows that modulo the extreme fragmentation regions, the enhanced hypron 
yield in SS and SAu can be understood in terms of the VENUS model double string 
hypothesis and in strong disagreement with extrapolated yields using HIJING from pp 
reactions. The factor of two enhancement of hyperon production is already needed in 
minimum bias pS in Fig. 3b, but is even more evident in Fig.5. The excess rapidity shift 
of the hyperons relative to non-strange baryons is also clear from Fig.5 since the peaks 
are shifted approximately one unit of rapidity further than those calculated for protons 
in Fig.4b. This enhanced hyperon transport is currently not explained by any of these 
models. In the preliminary Pb + Pb analysis|| a truly astonishing sharply peaked A 
distribution was suggested with dN(X)/dy = 23 ± 3. If conhrmed, A/p ~ 0.7 would be 
one of the most remarkable signature of novel phenomena in the heavy ion reaction. 

Since the MCM can account for at least the central region, non-strange baryon how 
up to SPS energies, it is of interest to extend the calculation up to RHIC energies, where 
the rapidity gap opens up to over 10 units. In Figure 6 we show the expected valence 
proton distribution expected at RHIC compared to present energies for Pb + Pb at b < 1 
fm central collisions. From this we estimate that the central baryon density that can be 
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conservatively expected at RHIC is 


3 dNP-P 1 
2 dy r p TiR 2 T 


~ Po (0.7 fm/r ) , 


(13) 


where p 0 ~ 0.17 fm 3 . 


5 Energy Degradation versus Baryon Number Trans¬ 
port 


A question related to baryon transport is the rapidity distribution of energy lost by the 
valence baryons. Do a few high rapidity pions carry away the missing baryon energy or 
is the energy shared between a large number of slower pions? Data on the forward veto 
calorimeter distribution shed useful light on this problem as was first emphasized in the 
analysis of WA80 data [[39(1 . 

In Figure 7 the veto calorimeter cross section is shown for NA35 S + S |h|] (Fig. 7b) 
and NA49 Pb + Pb (Fig. 7d ). The veto calorimeter measures energy in a narrow 
angular cone with A 6 ~ 0.3°. In both cases the well known horseback shape follows 
directly from Glauber geometry of spectator nucleons and is reproduced by VENUS and 
HIJING. However, the tail region at small Eyeto is sensitive to the energy degradation 
in central collisions. In Fig. 7 a,c the dependence of the contribution from collisions 
in the range b < 1 fm is illustrated in both VENUS and HIJING models. For central 
collisions, VENUS,in closer agreement with the data, displays greater energy degradation 
than HIJING. 

The correlation of the mean veto energy with the total multiplicity, participant nucleon 
number and impact parameter in both models is shown in Fig. 8. Fig. 8d demonstrated 
the equivalence of the Glauber multiple-collision geometry used in both models. However, 
the other parts show that the veto energy is systematically higher (and hence the energy 
degradation lower)in the HIJING model. Figure 9 shows the analogous correlations for 
the Pb+Pb reaction. 

While at first sight it may seem obvious that greater baryon stopping implies greater 
energy degradation, we show in Figure 10 that such a correlation cannot be taken a priori 
for granted given the uncertainty in the soft hadronization mechanism. In an older version 
of the diquark fragmentation scheme used in the version JETSET6.3|40| it was possible 
modify the baryon fragmentation through the parameter p = PAR( 52) in the LUDAT1 
common block. This parameter controls the momentum fraction of the junction J quark 
of a given diquark through 

f(xj) oc Xj( 1 - Xj) P , (xj) = — (14) 

3 + p 


Increase of the parameter p beyond the default value p = 1 obviously softens the diquark 
fragmentation function and leads to larger baryon rapidity shifts. Using the ATTILAf4l[ 


version of the FRITIOF model ||26|| it is possible to explore the consequences of varying 








this nonperturbative model parameter. The results are shown in Fig. 10. We note that 
HIJING1.3|28|] utilizes the new particle conventions of JETSET7.2p2|, which are incom¬ 
patible with the previous version 6.3. Unfortunately, in version 7.2 and higher, it is no 
longer possible to vary par(52). The once option (ihpr2(ll)=2 in HIJING) to switch to 
the JETSET “popcorn” fragmentation scheme unfortunately leads to even greater dis¬ 
crepancy with respect to pp —> pX data in Fig.l. 

It is seen from Fig. 10 that setting p = 10 for Pb + Pb reactions leads to a very large 
modification of the net proton distribution, from the default minimum to a maximum. In 
this case (xj) ~ 0.15. In comparison, the MCM model leads to (x) = (a/(l + a)) n ~ 1 /2 ~ 
(3/4) 3 /2 ps 0.2 for central Pb + Pb where the average number of collision per wounded 
nucleon is (z/) ~ 4.2. The MCM curve from Fig.4 is also shown for comparison. 

We find, therefore, that softening the diquark fragmentation can lead to similar baryon 
number transport as MCM and VENUS. Recall that VENUS incorporates also a diquark 
breakup component in its fragmentation scheme. The important point here is not that 
a variation of p can simulate more closely the baryon transport in VENUS and MCM, 
but that this model shows a priori no correlation between the energy degradation and 
baryon number transport can be assumed. In Fig. 10b the veto calorimeter cross section 
is essentially independent of p parameter. In this model the missing baryon energy is 
carried away into the veto calorimeter by a few high rapidity pions that fragment from 
the leading quark. The mechanism in VENUS that leads to greater energy degradation 
with increasing baryon rapidity shift is the dynamical assumption of the occurrence of a 
double rather than a single string in a fraction of the events that lead to large baryon 
rapidity shifts. From this analysis we conclude that the diquark mechanism of JETSET 
must be replaced in HIJING to enable the simultaneous reproduction of both the rapidity 
distributions and the veto calorimeter systematics. 


6 Conclusions 

One of our conclusions is that the nuclear stopping power, extracted from p+A data@ and 
extrapolated linearly via the MCM model to nuclear collisions, accounts quantitatively 
for the observed mid-rapidity baryon number transport up to central Pb + Pb reactions 
at 160 AGeV AGeV[l], [2, |3j]. The main difference between the valence baryon distribution 
functions in pA and BA arises from the variation of the Glauber geometrical probabilities 
of multiple interactions. The pion rapidity distributions are insensitive to the details of 
the energy degradation associated with the baryon number transport and are proportional 
to the wounded nucleon number. 

The forward veto calorimeter data is a sensitive probe of energy degradation and 
amplifies subtle variations of the angular distribution of fast pions. The HIJING model 
and its variants (with modihed diquark fragmentation via eq.([14|), cannot account for the 
veto cross sections. VENUS, on the other hand, reproduces both the baryon distribution 
and veto cross sections well. It remains to find out whether the novel double string 
mechanism assumed in VENUS is critical to both observations, or could the veto cross 
sections simply reflect an enhanced final state interaction effect in the spectator regions. 
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In any case, the strongly non-linear dependence of hyperon production definitely requires 
new dynamical mechanism that becomes operative already in p+S reactions. The VENUS 
model identifies this new mechanism with the double string formation. 

However, the large hyperon rapidity shift in the preliminary p + S and S + S data 
remain unaccounted for even with the present double string mechanism in VENUS 4.12. 
Not only is there a large enhancement of midrapidity hyperon production in p + S and 
S + S, but there is a rapidity mismatch between strange and non-strange baryons. Only 
the non-strange baryon transport is linear with A x / 3 within the present errors of the 
experiments. The provocative preliminary data of 0, |5j also point in the same direction. 

Hyperon production in Pb + Pb is however not entirely settled experimentally. The 
NA44 data in Fig. 4d include a substantial fraction of the net hyperon production while 
in NA49, the hyperon contribution has been subtracted out. The present experimental 
and systematic errors are too large to rule out a midrapidity maximum of the net baryon 
distribution. Our conclusion on linearity of baryon transport is thus subject to the above 
caveats. The present analysis in any case demonstrates the necessity of modifying the 
baryon fragmentation part of HIJING, FRITIOF, and DPM models. 
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Figure 1: Solid circles show the valence proton rapidity distribution (dn p /dy — drip/dy) in 
non-diffractive pp reactions at 400 GeV f38|, and open circles show the net proton distri¬ 


bution. The solid curve corresponds to the default Multi Chain model [||] fragmentation 
function eq. (§) with r p = 0.53. The dashed histogram is obtained using the ffIJINGl.3 
2?| code. 
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Figure 2: Baryon transport in minimum bias p,p-\-Ag reactions at 100 GeV. The solid and 
dashed histograms display HIJING and VENUS model results respectively as a function 
of Ay = ybeam — y with ybeam = 5.4. The smooth curve is the result of the Multi-Chain 
Model. The channel shown in part (b) isolates the target fragmentation region while (c) 
isolates the projectile fragmentation region. Part (d) shows the predicted A dependence 
of the proton rapidity distribution with a = 3 in the MCM. 
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Figure 3: Baryon number transport in minimum bias p + S, Au reactions]!] at 200 GeV. 
Curves are as in Fig.2. Here y is the laboratory rapidity with y > 0.5. In parts (a) and 
(c) the p — p distributions are shown while parts (b) and (d) correspond to the A — A 
distributions. Note that the data indicate significantly greater rapidity shifts hyperons 
relative to non-strange-baryons in contrast to both VENUS and HIJING calculations. 
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Figure 4: Comparison of central S + S at 200 AGeV (a,b) NA35|L], |20| and central Pb+Pb 
at 160 AGeV (c,d) (NA49|2|, |1^, |15[ (solid triangles), NA440 (solid squares)) data with 
calculations. Open circles and open squares are reflected data around mid-rapidity. New 
preliminary analysis^ of p — p are shown by solid dots in part (d). Solid and dashed 
histograms correspond to HIJING|2S[] and VENUS [2§| models, and continuous curves are 
predictions of the Multi-Chain Model f8|. 


16 






































rapidity 



"O 




Figure 5: Net Hyperon A — A rapidity distributions in central S + A at 200 AGeV (a,b) 
NA35[p]]. Solid and dashed histograms correspond to IIIJING|[28]| and VENUSQ models. 
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Figure 6: The energy dependence of the net positive baryon number transport in central 
Pb + Pb reactions based on the MCM. 
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Figure 7: The veto calorimeter cross section in S'S'[[T4|] and Pb-Pb [|15J collisions. Solid and 
dashed histograms refer to HIJING and VENUS models. 
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Figure 8: Dependence of the mean Eveto as a function of (a) the total multiplicity, 
(b) number of participant (or wounded) nucleons, and (c) impact parameter in S + S 
reactions comparing HIJING (solid) and VENUS (open) models. Part (d) shows the 
mean participant number as a function of impact parameter. 
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Figure 9: Dependence of the mean Eveto as a function of (a) the total multiplicity, 
(b) number of participant (or wounded) nucleons, and (c) impact parameter in Pb + Pb 
reactions comparing HIJING (solid) and VENUS (open) models. Part (d) shows the mean 
participant number as a function of impact parameter. 
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Figure 10: (a) Variation of baryon transport using the ATTILA [[4d 


version of Fritiof|26| . 

The effect of changing the parameter p = par( 52) in the jetset6.3ff(|] diquark fragmen¬ 
tation scheme is shown. The default (solid) with p — 1, and modified p = 5 dashed and 
p — 10 doted histograms are shown compared to the MCM prediction|| as in Fig. 4d. 
(b) The veto calorimeter cross sections and hence energy degradation are uncorrelated to 
the baryon number transport in this model. 
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